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We have investigated the motion of a single optically trapped colloidal particle close to a limiting wall at 
time scales where the inertia of the surrounding fluid plays a significant role. The velocity autocorrelation 
function exhibits a complex interplay due to the momentum relaxation of the particle, the vortex diffusion in 
the fluid, the obstruction of flow close to the interface, and the harmonic restoring forces due to the optical trap. 
We show that already a weak trapping force has a significant impact on the velocity autocorrelation function 
C{t) = {v{t)v{0)} at times where the hydrodynamic memory leads to an algebraic decay. The long-time 
behavior for the motion parallel and perpendicular to the wall is derived analytically and compared to numerical 
results. Then, we discuss the power spectral densities of the displacement and provide simple interpolation 
formulas. The theoretical predictions are finally compared to recent experimental observations. 

PACS numbers: 82.70.-y, 05.40.Jc, 87.80.Cc 



I. INTRODUCTION 

Understanding and controlling the motion of small col- 
loidal particles suspended in fluids and confined to small vol- 
umes is essential in many applications in microfluidics iQl] 
and biophysics |2]. Thermal fluctuations of the suiTound- 
ing fluid agitate the colloids giving rise to Brownian motion. 
Positioning the particle close to an interface or membrane 
and the direct optical observation of its trajectories allows 
in principle to use the bead's motion as a sensor for the 
chemical and physical properties of these two-dimensional 
surfaces. For example, the presence of surfactants lU, sur- 
face tension and elasticity 10, Et], or an adsorption layer |@] 
modifies the frequency-dependent mobility close to the in- 
terface. To establish Brownian motion as a local reporter 
of frequency-dependent surface properties a detailed under- 
standing of the complex interplay of the hydrodynamic flow 
and the colloidal thermal fluctuations is a prerequisite. The 
simplest effect known as surface confinement, namely the 
anisotropic reduction of the diffusion coefficient close to a 
bounding wall has been directly observed experimen- 

tally only very recently . To gain insight beyond the 
transport coefficients, one should record the time-dependent 
information, e.g. the mean-square displacement or the ve- 
locity autocorrelation, at times scales before the regime of 
simple diffusion is attained. Using weak optical trapping the 
quasi-free Browniari particle's trajectory can be monitored 
interferometrically 11 ill with a time resolution of a few mi- 
croseconds and a simultaneous positional sensitivity in the 
sub-nanometer regime lfl2i[l3ll . 

In the diffusive regime the momenta of the colloid are 
equilibrated with the surrounding fluid and are irrelevant for 
the Brownian motion. At very short times the particle's mo- 
tion is ballistic as expected by Newton's laws. Hence, only 
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after the colloid's initial momentum is transferred to the fluid 
the particle undergoes a random walk characterized by diffu- 
sion coefficients. The time scale where this transition occurs 
may be naively estimated by balancing the particle's inertia 
with the Stokes drag, Tp = mp/GTrjya, where nip denotes 
the mass of the colloid, a its radius, and jy the shear viscos- 
ity of the fluid. Since momentum is conserved it can only be 
transported by vortex diffusion giving rise to another charac- 
teristic time scale Tf = a^pf/i], where pf is the density of the 
fluid. The time a vortex emitted from the colloid to reach the 
wall separated at a distance h may be estimated by the char- 
acteristic time = h^pf/i]. The dynamical information on 
how the diffusive regime is reached, is encoded in the veloc- 
ity autocorrelation function (VACF) Cy (t) = (f y {t)v\^ (0)) 
for the motion parallel to the surface and similarly perpen- 
dicular to the wall. In bulk the VACF exhibits a power-law 
decay t^'^/^ at long times due to slow vortex diffusion. This 
long-time tail is also expected in the presence of a wall for 
times up to r„. For longer times this leading non-analytic 
behavior is canceled and a more rapid decay is expected, 
since the wall can carry away the momentum much faster. 
An algebraic decay t^^/^ for motion parallel and a t^^/^ 
long-time tail for motion perpendicular to the wall has been 
predicted by Gotoh and Kaneda lfl4ll extending an earlier 
work of Wakiya ifTsIl . However, the result for the perpendic- 
ular motion is erroneous, as pointed out by Felderhof lfl6ll . 
He succeeded also in calculating the frequency-dependent 
mobility for all frequencies in both directions and thus giv- 
ing an analytic solution for the motion close to the wall up 
to Fourier transform. Later, he generalized this work to the 
case of a compressible fluid lUTIl and to a second bound- 
ing wall O^. Computer simulations for a colloidal particle 
confined between two walls demonstrate that the long-time 
behavior of the VACF is strongly affected by the confine- 
ment CIQ. 

Recently, we have reported direct measurements of the ve- 
locity autocorrelation function of a single colloid immersed 
in water close to a bounding wall ll24ll . We have observed 
that the time-dependent VACF becomes anisotropic and ex- 
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hibits the non-algebraic tails that had been anticipated much 
earlier ifTiljl^ . To observe the particle for sufficiently long 
times close to a wall we used a trap constraining the Brow- 
nian motion even more. The optical trap introduces a har- 
monic restoring force — /-iTx, where x denotes the displace- 
ment from the trap center and K the spring constant. Ig- 
noring the wall a new time scale can be constructed = 
Gnrja/K characteristic of the positional equilibration in the 
trap. 

Here we supplement our experimental results with a 
derivation of the theoretical description of the Brownian mo- 
tion confined by a harmonic potential and a bounding wall. 
In particular, we discuss the complex interplay of the weak 
trap with hydrodynamic memory originating from the ob- 
structed vortex motion. We show that even if the trapping 
time Tk exceeds the characteristic time scales of particle mo- 
mentum relaxation Tp, fluid momentum diffusion Tf, and the 
wall-vortex reflection t„, by several orders of magnitude it 
still has a significant influence on the velocity autocorrela- 
tion function. Section |ll] provides a brief introduction into 
the theoretical framework of Felderhof, which is extended 
by a harmonic restoring force in Sec. |III1 An analytical 
discussion of the emerging long-time tails is presented in 
Sec. HV] followed by a numerical study of the VACF and 
the power spectral density of the displacement in Sec.[V] A 
comparison to our experimental results is shown in Sec. I VII 



II. FELDERHOF'S FRAMEWORK 

In this Section we review the theoretical basis elaborated 
by Felderhof iH for the motion of a single colloid in the 
vicinity of a wall. The fluid is treated as a continuum de- 
scribed by the Navier-Stokes equations, the colloid is mod- 
eled as an impenetrable sphere embedded in a viscous in- 
compressible fluid. At the particle's surface the usual no- 
slip boundary conditions are imposed. By the fluctuation- 
dissipation theorem the velocity autocorrelation function is 
related to the frequency-dependent admittance jzsll . hence 
it suffices to calculate the deterministic velocity response of 
the colloid to an external force. 

The external force R(t) causes the motion of both the par- 
ticle and surrounding fluid, which reacts by exerting a time- 
dependent drag force. Newton's second law for the acceler- 
ation of the sphere reads after a temporal Fourier transform 
— ia;(mp — mf)XJ^ = —F^ + R^^, where is the veloc- 
ity of the colloidal particle, iri j, its mass, and mf denotes the 
mass of the displaced fluid. The total induced force F^j cor- 
responds to the frequency-dependent drag force on the par- 
ticle up to the acceleration force — iwmfUt^ of a rigid sphere 
of fluid of equal radius. We shall employ the point particle 
limit where inhomogeneities of the flow on the scale of the 
colloid are ignored. 

To calculate the total induced force one maps the problem 
in the presence of a bounding surface to a corresponding one 
in infinite space and employs the generalized Faxen theo- 
rem of Mazur and Bedeaux 12611 : For an unbounded fluid the 



frequency-dependent induced force F^ is 
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with ({uj) ~ 67r7]a(l + \/— icjrf) and denoting the unper- 
turbed flow evaluated at the position of the particle Tq. The 
branch cut is chosen such that \/—iujT{ = (1 — i)y^ujT{/2. 
Then consider as acting flow 
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where Vi^(r), vocj(r) denote the solutions of the Navier 
Stokes equation in response to a point force acting at 
ro in the presence of a bounding surface and in the infinite 
space. Clearly, v^(r) satisfies the homogenous Stokes equa- 
tions. The key idea is that Vtj(r) can be interpreted as a flow 
in infinite space resulting from the force R^^ with (r) as 
an externally acting flow. Yet, for flows in infinite space the 
generalized Faxen theorem applies and the induced force can 
be calculated easily. 

By linearity of the Navier-Stokes equation, the fluid re- 
sponse to a point force R^ acting at Vq is obtained by 

v^(r) = G'(r,ro) • Rc^ , voi^(r) = Go(r - ro) • R<^ , 

(3) 

where G(r, rp) and G'o(r — Tq) are the corresponding (ten- 
sor) Green's functions. Then the acting flow at the position 
of the particle reads = F{ro, uj)Il^, where 



F{ro,Lj) = lim (G(r,ro) 

r->ro 



Go(r,ro)) , 



(4) 



defines the reaction field tensor. Its frequency dependence 
only enters via the ratio h/6 = u]Tf/2, where h is the 
distance to the wall and 6 ~ yj2ii/ pfuj denotes the skin pen- 
etration depth. Combining Newton's second law, Faxen's 
theorem, and the reaction field approach one finds for the 
force balance 



[U^-J^(ro,w)R„)]+R„.(5) 



Solving this equation for the particle velocity yields U,^ = 
3^(w)Rtj, where the response function y{'jj) is called the ad- 
mittance tensor and corresponds to a frequency-dependent 
mobility. 

Consider first the case where the wall is infinitely far away 
and the reaction field tensor vanishes. Then one finds for the 
admittance for infinite space 



3^0 (w) = [-iojm* + C(w)]' 



(6) 



Here to* = + mf/2 can be interpreted as the effective 
mass of the particle, since for a particle moving in an ideal 
fluid with constant velocity U the total kinetic energy in- 
cluding the dragged fluid is given by to*U^/2 i27ll . 
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The general result in presence of a bounding wall can then 
be expressed as 



close to a confining wall with the externally acting harmonic 
restoring force. 

Following the chain of arguments of Sec. we include 
the harmonic restoring force in Newton's second law 



1 + GTTjya ( 1 + 1/— icjTf + 
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F{ro,io) 



(7) 



Thus, all modifications due to the bounding wall are de- 
scribed by the reaction field tensor By symmetry only the 
motion parallel and perpendicular to the wall 3^||(ro,aj), 
y±{ro, Lu) have non-zero components. Felderhof succeeded 
in providing a full analytical result for the frequency de- 
pendence of the corresponding reaction field tensors. We 
have verified the sophisticated results for the parallel mo- 
tion, Eq. (3.9) in lUol and the perpendicular motion I28i1 . 
For further analysis, in particular, to evaluate the long-time 
tails in the velocity autocorrelation function, it is sufficient 
to know their corresponding low-frequency expansions 
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where v ~ (— iajT„)^/^. Specializing to the stationary case, 
cj = 0, and inserting into Eq. O, one easily obtains the 
zero-frequency admittance 3^|| 0) recovering Lorentz's 
result 0] for the mobility close to a wall 



H = MO 
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where fig — ^/ GT^rja denotes the mobility in bulk. 



III. INFLUENCE OF A HARMONIC TRAP 

In the framework presented above it is implicitly assumed 
that the distance h between the particle and its bounding wall 
is time-independent. To realize this situation experimentally 
and to prevent the particle from leaving the detection re- 
gion a trap is needed, which should be included to complete 
the model. Already in bulk, the motion of a harmonically 
bound Brownian colloid differs drastically from a free parti- 
cle, since at long enough times the restoring forces limit the 
particle's mean-square displacement to a finite value. The 
very first attempt by Uhlenbeck and Ornstein ll29tl to model 
the interplay of inertia of a particle. Stokes friction and har- 
monic restoring forces dates back to 1930 but ignores hy- 
drodynamic memory effects. Only much later, Clercx and 
Schram Isoll incorporated the fluid inertia generalizing the 
VACF for the free motion. In particular, they found that 
the long-time anomaly changes from a t~^^^ behavior for a 
free particle to a much more rapid decay according to 
Here, we combine the ideas of Felderhof for the free motion 



ia;(mp — mf)\Ji^ 



(10) 



Here K denotes the spring constant of the optical trap, which 
we assume to be isotropic. The displacement of the particle 
relative to the trap center at ro will be eliminated in favor of 
the particle's velocity via = —iujx^. Since Faxen's law 
and the reaction field tensor are unaffected by the presence 
of the trap, one can again eliminate the total induced force 
and merely supplement the force balance, Eq. by the 
harmonic force 



K 

-iuj{mp - mf)lJ^ - — U, 



C(w) - -i ujirif 



In terms of the admittance tensor the previous relation may 
be written as 



— U„ + R, 

lUJ 



(11) 



with the following interpretation: The particle's velocity is 
still determined by the admittance tensor without trap, pro- 
vided now the net force is considered as driving the system. 
Solving for the velocity = J^*^'^^ (ro, a;)R(^, one finds 
that the harmonic potential modifies the admittance accord- 
ing to the simple rule 



(12) 



This result includes Felderhof's result for zero trapping, 
K = 0, SiS well as Clercx and Schram Isoll if the wall is 
infinitely far away, 3^(ro, w) — >■ 3^o('^)- Again by symmetry, 
the matrix y'^'''> (rg, uj) is diagonal with elements J^i^i*^-* (rp, co) 

and J^j'"' (ro , w) for the force R-^ parallel or perpendicular to 
the wall. 

The trap introduces a new characteristic time scale into 
the problem: balancing the harmonic restoring force Kx 
with the zero-frequency Stokes drag Girrjav, one obtains 
Tk = 6m]a/K = I/iiqK. The reference drag force has been 
chosen as the one acting in bulk, although it is clear that the 
presence of a boundary suppresses the hydrodynamic fric- 
tion. The physical interpretation of the dependences is that 
a stronger trap will lead to a faster relaxation to equilibrium 
resulting in a reduction of Tk, whereas an increase of the 
Stokes drag induced by a larger viscosity leads to slowing 
down of the equilibration process. 

In the derivation we have assumed that the trapping is 
isotropic, i.e. the spring constants are identical for all direc- 
tions of the displacement. In case the trapping potential be- 
comes anisotropic, the restoring force is still given by Kx^ 
where the spring constant K has to be interpreted as a sym- 
metric matrix. Then Eq. ( fT2] i still holds, and if the principal 
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axes of the trap includes the direction perpendicular to the 
wall, the matrix inversion is achieved by inverting the diag- 
onal elements. 

Let us compare Eq. ( fT2] i to the well known high- and low- 
frequency limits. For high frequencies the trap becomes in- 
creasingly irrelevant and the response is dominated by the 
inertia of the particle and the displaced fluid. In the low- 
frequency regime the harmonic potential suppresses the ad- 
mittance reflecting the fact that no static external force can 
induce a stationary motion of the particle in confinement. 



IV. ANALYTICAL DISCUSSION OF THE VELOCITY 
AUTOCORRELATION FUNCTION 

Here we focus on the velocity autocorrelation function 
(VACF) Ci\{t) = ('U||(t)?;||(0)) for the motion parallel and 
perpendicular to the wall, C±{t) = {v±{t)v±{0)) . The 
VACF are connected to the admittances of the previous Sec- 
tion via the fluctuation-dissipation theorem 
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where ksT denotes the thermal energy. This relation holds 
for all cases, i.e. parallel and perpendicular to the wall, with 
and without trap. To simplify notation we have suppressed 
the dependence on the position rp. 

Due to trapping the particle samples only a region close to 
the center of the trap and according to a stationary ensemble 
given by the Gibbs-Boltzmann measure. In particular, time- 
dependent position correlation functions like (x(t) • x(0)) 
are well defined. In the derivation for the admittances, the 
position of the particle was assumed to be fixed at Tq. How- 
ever, since the particle fluctuates and takes excursions from 
the center of the trap of typical magnitude ^Jk^TJK, the 
spring constant of the trap has to be strong enough in order 
to render these fluctuations negligible compared to the dis- 
tance to the wall h. 

The short-time evolution of the velocity autocorrelation 
is inferred from the high-frequency behavior of the corre- 
sponding admittances. Obviously, the harmonic restoring 
forces do not affect the initial value and they are still given 
by 
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Without a bounding surface C^^\t = 0) = kBT/m* |13C 
is recovered. The equipartition theorem suggests that the 
initial value should read ksT/rjip for all cases. However, 
as discussed already by Zwanzig and Bixon llzsll . there is an 
additional rapid decrease of the VACF if the finite compress- 
ibility of the fluid is taken into account. 

The low-frequency expansion of the admittance exhibits 
non-integer powers, which by tauberian theorems llsill cor- 
respond to algebraic long-time decays in the VACF. The 
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TABLE I: Leading non-analytic term in the low-frequency expan- 
sion of the dimensionless admittances &TiT]ay{ijj) for the free case 
(0), perpendicular (_L) and parallel (||) and again in the presence of 
the trap (k, k_L, k||) 



leading long-time behavior of the normalized velocity au- 
tocorrelation, C{t — > oo)/C(0) ~ (r/t)", can be inferred 
from the leading non-analytic contribution in the admittance 
Y{u} -> 0) ~ C(0)T(-iwT)"-ir(l - aj/ksT. The phys- 
ical origin of the long-living correlations lies in the long- 
range fluid flow that is generated by momentum conserva- 
tion in the Navier-Stokes equation. At low Reynolds num- 
ber, transverse momentum can only be transported away by 
vortex diffusion, hence there is a time-dependent growing 
length scale R{t) ^ {rjt/ pi)^/"^ characterizing how far mo- 
mentum has penetrated into the fluid. In the simple case of 
unconfined free motion, the particle's initial momentum is 
then shared with a fluid volume ^ R{tY leading to a long- 
time behavior of C(t)/C(0) ^ t"'^/^. For the free motion 
close to a bounding surface, the leading term ~ ^""^/^ is 
canceled and a more rapid algebraic decay t~^/^ lfl6ll is ex- 
pected. Interestingly, the same algebraic decay occurs also 
in the disordered Lorentz gas Il32l - B4ll . a simple model for 
transport in porous media. There the long-time memory 
arises since the particle remembers the presence of an ob- 
structing wall for arbitrarily long time. Confining the motion 
of a colloidal particle in bulk by a harmonic potential leads 
to a more rapid decorrelation according to i~ ""/^ Isoll . In the 
case of a trapped particle close to a wall, a series expansion 
3^11 (w), J^^*^'' (w) in powers of the frequency reveals a lead- 
ing non-analytic term io^/"^ resulting in a ~' 
obtain this result it is sufficient to know the reaction field 
tensors (tq, w), F±{tq^lu) including 0{lli) terms. Table|l] 
summarizes the leading non-analytic behavior of the admit- 
tances including all prefactors. The corresponding long-time 
tails in the normalized velocity autocorrelation function are 
displayed in Table [III where for the walls only the leading 
term in T„/Tf is shown. 

For the motion perpendicular to the wall, the true long- 
time behavior is masked by the next-to-leading term. This 
has been pointed out by Felderhof lfl6ll for the trap-free mo- 
tion and here we supplement the analysis also for the trapped 
dynamics. In order to derive the next-to-leading term in the 
low-frequency expansion the reaction field tensor F±_ (w, rg) 
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TABLE II: The leading long-time behavior of the normalized ve- 
locity autocorrelation function C{t ^ oo)/(7(0) for the free case 
(0), perpendicular {IS) and parallel (||) and again in the presence of 
the trap (k, k_L, k||). For the cases of abounding wall, only the lead- 
ing contribution in r„/rf is displayed. The mass ratio is encoded in 
the constant B = (9rp/rf + l)/180r. 



has to be evaluated including the order (— iwrf)^/^. Then 
the two leading nonanalytic terms in the far-field expansion 
Tw/Tf > 1 read 
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without harmonic restoring force, and 
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including the trap. The second term becomes negligible with 
respect to the first one only at frequencies below uj ^ Tf/ t^. 
For the corresponding time correlation functions this implies 
a long-time behavior according to 
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The amplitude of the leading term depends on the ratio of 
the mass of the particle to the displaced mass of the fluid 
or equivalently on the ratio Tp/rf = 2mp/9TOf. For typical 
experimental conditions, e.g. a silica sphere in water, the 
amplitude is negative and one should expect the correlation 
function to approach zero from below. The presence of the 
factor Tw/ Tf renders the amplitude of the subleading term to 



a large quantity. Balancing both terms yields an estimate 
/Tf for the time one has to wait in order to observe the 
true long-time behavior Hence the perpendicular motion 
without trap exhibits a zero crossing on a time scale much 
longer than the naive guess r„. We shall see in the next Sec- 
tion that even a weak trap has strong implications for this 
intermediate time behavior 



V. NUMERICAL RESULTS 

A. Velocity autocorrelation function 

The inverse Fourier transformation cannot be calculated 
analytically, since the admittance tensors are not simple ele- 
mentary functions. A numerical Fourier transformation can 
easily be performed, and since the time-correlation functions 
are real and even in time, a Fourier-Cosine transformation is 
sufficient 



C{t) 



2kBT 



dcj cos(a;ORe [y{u)] . (19) 



Since the VACF exhibits long-time power-law behavior the 
admittance has to be sampled over many decades in fre- 
quency. A conventional Fast Fourier Transform (FFT) is 
rather ill-suited to achieve this and would also calculate 
C{t) on an equidistant time grid not convenient for our 
purposes. Hence, we apply a simple modified Filon algo- 
rithm ifssll with typically 10^ frequencies per decade cov- 
ering 30 decades to calculate C{t) on a logarithmically 
equidistant grid containing 100 data points per decade. We 
employed Mathematica® to evaluate the various special 
functions appearing in the reaction field tensors and spliced 
them together with a high-order series expansion at low and 
high frequencies. We have checked that our numerical data 
reproduce the predicted tails over several orders of mag- 
nitude, with the exception of c'"^\t), where the next-to- 
leading tail dominates the figures. 

The velocity autocorrelation for the motion parallel to a 
wall is displayed in Fig. [TH for a colloid gradually approach- 
ing the bounding surface. Without trap the correlation func- 
tions remain positive for all times. The long-time behavior 
exhibits the t^^/^ decay with an amplitude that depends sen- 
sitively on the distance to the wall. For the weak harmonic 
restoring force under consideration here, the VACF follows 
the free behavior up to intermediate times; at longer times 
the VACF is strongly influenced by the trap. In particular, 
the VACF exhibits two zeros as is the case also without the 
bounding wall. Note that the first zero is by orders of mag- 
nitude earlier than the characteristic time scale of the trap 
Tk, however much later than the Langevin momentum relax- 
ation time Tp. Between the two zeros a negative flat plateau 
emerges that extends to time scales longer than Tk. The very 
late decay is governed by the tail of Table Ull however 
its onset exceeds the range of the figure. 

The VACF in the direction perpendicular to the wall is ex- 
hibited in Fig. [TJ). For a typical experimental particle with 
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FIG. 1: (Color online) Double-logarithmic representation of the 
normalized VACF for the motion parallel (a) and perpendicular 
(b) to the wall. The parameter Tp/rf = 0.5 corresponds approx- 
imately to a silica bead in water. Full lines correspond to weakly 
trapped particles rf/rk = 10^'', whereas the dotted lines are with- 
out trapping (rf/rk — 0). The wall is gradually approached, 
Tw/tt — 16, 8, 4, 2, corresponding lo h/a = 4, 2.83, 2, 1.41, and 
the initial decay shifts to the left. The straight lines are guides to 
the eye for the t^''^^ power law. 



momentum relaxation time Tp = 0.5rf and without confin- 
ing potential the curves exhibit a single zero. In this case 
the amplitude of the leading long-time behavior be- 
comes negative and the next-to-leading order i"^/-^ domi- 
nates at an intermediate time interval, see Eq. ( fTTj i. The cor- 
relation functions without trapping thus approach zero from 
the negative side by an algebraic decay that is fairly insen- 
sitive to the distance from the wall. The zero on the other 
hand shifts rapidly to shorter times as the wall is approached. 
From Eq. ( [TtI i one infers an asymptotic scaling behavior 
~ T^/Tf for the zero far away from the wall. For the two 
curves corresponding to the distant particle, Tw/rf = 8, 16, 
the tr'^^^ decay is visible for intermediate times close to the 
zero, as has been reported earlier lfl6ll . For a weak trap 
Tf/Tk = 10^^ the curves follow the ones of the uncon- 
strained motion down to a signal of 10^^ and start to de- 
viate strongly at later times. Again a flat negative plateau 
characteristic of the weak restoring force is attained. The 
decay from this intermediate plateau occurs at times large 



compared to followed by a rapid algebraic decay t~^/^ 
without passing through another zero (not shown). For dif- 
ferent parameters r^, t„ one finds also a scenario with three 
zero crossing. 

To study the hydrodynamic memory effects for the col- 
loidal motion close to the wall, we choose a trapping po- 
tential that is as weak as possible. Then the time scale 
Tk = Gnrja/k is much larger than the remaining time scales 
of the problem Tp , Tf, and r„ and one is tempted to argue 
that the trap is irrelevant in the regime of interest. How- 
ever, as is discussed in Fig. [T] the trap has significant im- 
pact on the velocity correlation functions even if the char- 
acteristic time scales differ by orders of magnitude. The 
reason for such a behavior is twofold: First, the hydrody- 
namic memory leads to a scale-free power-law long-time de- 
cay and the parameters Tp, Tf, and determine merely the 
amplitude of the algebraic behavior rather than a character- 
istic decay time. Second, for the point particle limit to apply 
accurately, h/a > 3, and the VACFs become small at times 
t > where the tail t^^/^ is expected to set in. Thus al- 
ready a weak trap has a strong influence on the signal where 
the interesting hydrodynamic memory effect dominates the 
VACFs. 



B. Power spectral density 

An alternative way to investigate the interplay of fluid 
inertia and a trap is to focus on the power spectral den- 
sity, ail approach that has been pursued by Berg-S0rensen 
et al. 113a, 13711 . There the fluctuating position x(t) of the 
bead in a finite time interval [— T/2,T/2] is decomposed 



into Fourier modes xt{oj) = ^(^) S:xp{iujt)dt where 

the angular frequencies ut are integer multiples of 27r/T. 
For long observation times, T — > oo the power spectrum 
S{uj) — (^\xj{Ld)\^'^ /T becomes a quasi-continuous func- 
tion for frequencies oj ^ 2tt/T. Since the corresponding 
velocities fulfill vj{lo) — —\ujxj{ljj) up to irrelevant bound- 
ary terms, the power spectrum can be also expressed as 
ijj'^S{ijj) = (|wt(w)P) /T. Thenwiflifliehelpofflie Wiener- 
Khinchin theorem aj^5'(aj) = /(u(t)u(O)) cxp(ia;t)dt and 
the fluctuation-dissipation theorem Eq. ( fT3] l. the power spec- 
tral density can be obtained from the admittance 



S{u) 



2kBTRe [y'-'^Huj)] 



(20) 



Since the trap modifies the admittance according to Eq. ( fT2] l. 
we find for the power spectrum 



2kBTRe [y{u 



(uj Re [y{uj)-^y + {cj\m[y{u 



KY 



(21) 

This result is again valid for all cases, i.e. in bulk, paral- 
lel and perpendicular to the bounding wall. For the zero- 
frequency limit in the case of a limiting wall we find 



^ll(O) 



2kBT 



S±{0) = 



2kBT 



(22) 



7 



using the zero-frequency limit of the admittance at the wall. 
Without trap, the power spectrum S{uj) would diverge for 
a; —> as can be seen from Eq. (|20| |. reflecting the fact that 
the particle can take excursions without bounds. The 
dependence on the trap strength arises from two aspects. 
First, the motion of the particle is harmonically confined 
leading to equilibrium fluctuations (x^) = ksT/K. Sec- 
ond, the time scale to reach equilibrium may be estimated by 
balancing the restoring force yielding l/piK. For the bulk 
motion this is identified with = Q-K-qa/K = 
whereas close to the wall the reduction of the mobilities 
, iJ,± due to the obstruction of the hydrodynamic flow has 
to be taken into account. 
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the particle, the power spectra assume a Lorentzian shape 

2kBT 1 



(23) 



where ^ = /io,/i||,/i^ is the corresponding mobility. For 
the weak traps employed here, this gives an accurate repre- 
sentation of the power spectral densities for dimensionless 
frequencies of order unity. 
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FIG. 3: (Color online) Power spectral density with respect to a 
Lorentzian, Eq. l l23b . in the frequency regime where fluid inertia 
plays a role for the parallel (a) and perpendicular motion (b). The 
parameters correspond to the ones of Fig. [2] Also included is the 
harmonic oscillator prediction according to Ornstein and Uhlen- 
beck. 



FIG. 2: (Color online) Semi-logarithmic representation of the 
power spectral density for a weakly trapped Brownian particle 
Tf/rk = 10"'' close to a limiting wall for the ratio Tp/rf — 0.5 
((a) lateral direction and (b) vertical direction). The distances to 
the wall are the same as in Fig. [T] The normalization is cho- 
sen as the zero-frequency limit of the power spectrum in bulk 
Soii^ = 0) = 2kBT/fioK'^. The increase of the initial value is 
due the wall according to Eq. l i22t . 



At the scale of the trap relaxation rate 1 /r^, the main fea- 
ture of the power spectral density as displayed in Fig |2] is 
the saturation at a height given by Eq. ( l22l i at low frequen- 
cies. A rapid decrease of the power spectrum is manifest at 
higher frequencies. Ignoring inertial effects of the fluid and 



To highlight the deviations from the simple Lorentzian, 
we display the ratio S{uj) / S^^^^ {ui) for the parallel and per- 
pendicular motion in Fig. [3] The differences arise at the fre- 
quency scale 1/Tf, where the motion of the hydrodynamic 
vortex sets in. At higher frequencies the full power spec- 
trum decays more rapidly than a Lorentzian due the iner- 
tia of the particle and the fluid. Interestingly, the deviations 
are strongest in bulk and they fade out only slowly as the 
frequency is decreased. The wall suppresses these inertial 
effects consistent with the notion that the rigid interface car- 
ries away part of the particle's initial momentum. Hence 
the obstruction of the vortex pattern, manifest in the time- 
dependent VACF as reduction of a to t^^/^-tm\, is also 
visible in the power spectral density by a faster approach to 
the Lorentzian shape for low frequencies ujTv, < 1. Note that 
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for the perpendicular motion the power spectral density be- 
comes larger than a Lorentzian very close to the wall. This 
phenomenon is related to the change of sign of the long-time 
anomaly of the t^^/"^ in Eq. (fTTT i. The overshoot disappears 
if denser colloidal particles are used, i.e. Tp/rf is increased. 
However, the vanishing of the overshoot does not coincide 
with the sign change of the power-law tail at Tp/rf = 5/9 
but occurs at larger ratios of approximately 1.5 for the pa- 
rameters used here. For comparison we have also included 
in both panels the harmonic oscillator prediction for the bulk 
motion following Ornstein-Uhlenbeck 

^'^> A^oi^Ml-wVpr02 + (wrk)2' ^^^^ 

neglecting both, the added mass due to the displaced fluid, 
as well as the vortex motion. As can be inferred from Fig.|3] 
the wall suppresses the hydrodynamic memory, and for the 
parallel motion the power spectral densities are typically in 
between the bulk behavior where hydrodynamics is included 
and the Ornstein-Uhlenbeck shape. 

For the weak trapping regime under consideration, the 
time scale is much larger than the other ones. Then one 
may assume that the processes of the relaxation of the mo- 
mentum to the surrounding fluid and the relaxation of the 
position in the trap are decoupled, which suggests to approx- 
imate the power spectral density as 

5H«5(^)m5^M, (25) 

where again the appropriate mobilities /i = /ip, M|| , and 
similarly for the admittances y{uj) have to be inserted for 
the different cases. We have checked that for the param- 
eters used here this constitutes an excellent approximation 
for all frequencies. Consequently, Fig.[3]essentially displays 
the real part of the admittances for different distances to the 
wall. The reason that the curves superimpose is that the zero- 
frequency limit serves as a background that overlaps the in- 
teresting non-analytic behavior discussed in Sec.|IV] Con- 
sequently one can not use Eq. (IZSl l as input in the numeri- 
cal Fourier transform since it does not result in the correct 
long-time behavior Concluding, although the power spec- 
tral density allows for simple approximations on different 
frequency scales, their deviations encode subtle correlations 
that are manifest and more directly accessible in the time- 
dependent VACF. 

VI. COMPARISON TO EXPERIMENTS 

Recently, we have performed experimental tests for a 
small silica sphere trapped by a laser focus in the vicin- 
ity of a surface ll24ll . The ttajectory of the particle is 
measured interferometricallv lUlll with a spatial resolution 
in the subnanometer range II12I1 . We employ an infrared 
(A = 1064 nm) diode-pumped Nd:YAG laser (IRCL-500- 
1064-S, CrystaLaser, USA) with a cw output power of 500 
mW. The beam is first expanded 20 x and then focused using 



a)| 



b)| 




FIG. 4: (Color online) (a) 3D lateral view of the experiment; 
spheres are drawn to scale. A silica particle of radius a =1.5|im 
trapped by the laser focus is placed next to the surface of a sig- 
nificantly larger silica sphere. This 100 [jm sphere is immobilized 
between the two coverglass surfaces of the sample chamber (b) 
Optical image of the probing particle's position relative to the wall 
created by the big sphere. The 3 \im probing particle was placed at 
a distance h =1 1.5 |im away from the 100 |jm sphere's surface and 
gradually approached. The velocity correlation functions, as well 
as the diffusion coefficients for the motion parallel and perpendic- 
ular to the wall are measured. 



a 63 X water-immersion objective lens of a numerical aper- 
ture NA=1.2 resulting in a stable optical trap for the colloid. 
An InGaAs quadrant photodiode (G6849, Hamamatsu Pho- 
tonics, Japan) is placed in the back focal plane of the con- 
denser lens recording the modulation of the optical power 
due to the displacement of the particle near the beam fo- 
cus. The photodiode signal is amplified and digitized using 
a data acquisition card with a dynamic range of 12 bits. The 
detected positions contain N = 10^ points separated by 2[rs, 
which corresponds to a sampling rate of fs = 500 kHz and 
a recording time of tg ~ 20 s. 

The Brownian particle is a silica sphere of radius a ~ 
1.5 |am and mass density pp — 1.96 g/cm'' immersed in wa- 
ter (pf = Ig/cm'^, 77 = 10~'^Pa-s) at ambient temperature. 
The characteristic time scale corresponding to the particle's 
inertia Tp = mp/67rr/a = 2a^pp/9?7 = 1 |is is about half the 
one of the fluid inertia Tf = a^pf/77 = 2.25 ps. Using optical 
tweezers, we approach the particle incrementally towards a 
sphere of diameter lOOpm, much larger than the size of our 
particle. Then sufficiently close to the large sphere the cur- 
vature can be ignored, and the particle undergoes Brownian 
motion close to a planar wall around the trap's center, see 
Fig-S] The large sphere is immobilized since it is in close 
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FIG. 5: (Color online) Log-log plot of both normalized VACF, 
C||(*)/C||(0) and Cx(f)/Cx(0) for a sphere (a =1.5 |jm, Tp = 
1 |js, Tf = 2.25 |js) trapped in a weak optical potential {k ~2 |jN/m, 
Tk = 14 ms). The increasingly anisotropic VACF is measured at 
three distances from the wall {h= 9.8, 6.8 and 4.8 |im, correspond- 
ing to T„ = 96, 46, 23 |is, respectively). Positive correlations are 
represented by full symbols, negative ones by open symbols. The 
characteristic power-laws are represented by thick lines as guide to 
the eye. 



contact with the two coverslides of our fluid chamber (size 
w 2 cm X 0.5 cm and thickness « 100 pm). In this set-up we 
have equal sensitivity for the motion parallel and perpendic- 
ular to the surface under consideration, since both directions 
are perpendicular to the axis of the laser beam. The sample 



is mounted onto a piezo-stage, and the 100 [im sphere can be 
positioned at a distance h relative to the trapped particle by 
moving the piezo-stage in all three dimensions with a preci- 
sion of « 1 nm. 

To study the hydrodynamic memory effects induced by 
the presence of a wall rather than the motion due to the trap 
confinement, we have applied the weakest trapping force 
possible. Yet, the optical trap should still be strong enough 
not to loose the particle from the laser focus during the ex- 
periment. Hence, we have optimized the trapping strength in 
order to suppress the effects of the trap on the VACF without 
loosing too much sensitivity at the detector A series of ex- 
periments revealed that k w 2fiN/m fulfills all these require- 
ments. Note that the corresponding time scale = 14 ms 
exceeds the parameters of the bulk motion Tp, Tf by four or- 
ders of magnitude and the one of the wall by more than 
two orders of magnitude even for the farthest distance stud- 
ied. However, this weak trapping force has still significant 
influence on the VACF at the time scale where the hydrody- 
namic memory becomes apparent. 

We have recorded the particle's position r|| (t„), rx(i„) 
at equidistant instants of time tn '■= nAt, At := l/.fs, 
n e No, and derive coarse-grained velocities as := 
[r\\{tn + At) — r||(i„)]/At. The velocity autocorrela- 
tion functions are then evaluated as time-moving averages 
C\i{t) = N-^J2n=i^Min + t)Aw||(i„), and similarly 
for the motion perpendicular to the wall. The normalized 
VACFs C||(t)/q|(0) and C±{t)/C_L{0) are displayed in 
Fig.|5]for different distances h from the wall. For large sep- 
aration h ~ 9.8 pm both VACFs coincide within our error 
bars down to where the signal reaches the level of 1%. This 
regime is clearly dominated by bulk behavior and, in par- 
ticular, the well-known algebraic decay t^^^^ characteristic 
for the unconstrained vortex diffusion is recovered. At the 
time scale = Pfh? /rj the vortex generated by the ther- 
mal fluctuations of the Brownian particle reaches the wall 
and the correlation functions split as the motion becomes 
anisotropic. The wall leads to a more rapid decay and for 
the parallel motion a power law t^^/^ enters the observa- 
tion window . The signal of the perpendicular motion is an- 
ticorrelated at these time scales. Later times lead to only 
weak signals that we cannot resolve within our noise level 
of 10^**. For the VACFs the trap manifests itself in two ze- 
ros in the parallel motion, the second of which is outside of 
our observation window. The perpendicular motion exhibits 
a first zero which is shifted to earlier times by the presence 
of the optical trap — a second zero induced by the harmonic 
restoring forces remains unobservable. Since for silica in 
water Tp/rf = 0.44 < 5/9, the theory ignoring the trap also 
predicts an anticorrelated signal in the long-time behavior 
for the perpendicular motion. However, the negative signal 
we observe is dominated by the restoring force of the optical 
trap. Approaching the wall to h = 6.8 |jm and 4.8 |jm, the 
characteristic time r„ = 46 |js, respectively 23 ps, decreases 
and the splitting into parallel and perpendicular motion shifts 
to earlier times. For the parameters chosen, the bulk behav- 
ior is dominated by the t~'^/^ tail, and the splitting sets in at 
higher values of the correlation functions rendering it easier 
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to observe. The crossover in the parallel VACF from the bulk 
dominated behavior ^^■^/^ to the wall dominated algebraic 
decay t^^/'^ becomes more and more pronounced. The zero 
due to the trap is dragged to shorter times since the surface 
constrains the colloid's motion increasingly. The perpendic- 
ular component C^{t) decreases more rapidly and practi- 
cally enters our noise floor at the expected zero induced by 
the trap. The asymptotic expansion for the long-time behav- 
ior of (t) suggests that a power law of with positive 
amplitude should be present at intermediate times. Such a 
behavior may be inferred from the data for a narrow regime 
of times, but the presence of the trap suppresses its ampli- 
tude by a factor of two. 

Next, we define the time-dependent diffusion coefficients 

D\\{t)^ f C\\{t')dt' , (26) 

and similarly for the motion perpendicular to the wall. From 

the recorded time series, D\\ {t) = Y^n=i [''II + ^ 
r|| (t„)]?j|[ (i„) is directly evaluated. We have checked that 
this gives the same result as integrating the VACF. The dif- 
fusion coefficients D\\{t),D±{t) corresponding to the dis- 
tances h = 9.8, 6.8 and 4.8 pm of Fig.|5]as well as the mo- 
tion in bulk {h ~ 37.8 |jm) are displayed in Fig. |6] The 
most prominent feature is a plateau extending to the time 
where the trapping becomes effective. The height of each 
plateau approaches the diffusion coefficients parallel and 
perpendicular to the wall that are obtained by the Einstein- 
Smoluchowski relation from the mobilities, D\\ = ksTji^^, 
D± = ksTj^L^. The reduction of the zero-frequency mobili- 
ties due to the wall in Eq. (|9|l ignores higher order correction 
terms m a/h which can be calculated exactly i^. For the 
distances studied here h/a > 3 the point particle limit is 
accurate within 2% and the correction terms can be safely 
ignored. Since the VACF exhibits long-time tails, the ap- 
proach to the plateau is slow. Ignoring the trap, we expect 
asymptotically for i — > oo in leading order in a/h 

where D ~ fc^T/GTrr/a denotes the bulk diffusion constant. 
Even for the weak trapping used in the experiment the max- 
imum in the time-dependent diffusion coefficient deviates 
from the asymptotic value by a few percent. The maximum 
corresponds to the zero crossing in the VACF, which for the 
parallel motion is entirely due to the trap. At longer times the 
diffusion coefficient (t) decreases and eventually reaches 
zero. The diffusion constant is given by the Green-Kubo 
relation _D|| ~ C\\{t')dt' or the long-time limit of the 
time-dependent diffusion coefficient. Hence, at long times 
no diffusion occurs, which is consistent with the observation 
that the harmonic restoring forces localize the particle at suf- 
ficiently long times. Then the data at long times are sensitive 
to the precise value of the trap time r^, and we have used this 
observation to optimize the fit for the trap stiffness. Since 
the laser power is held at a constant value for all distances. 



the trapping time should be identical for all data. Figure |5] 
shows that the terminal decay to zero can indeed be fitted by 
a single t^, however for the closest distance h = 4.8 [im the 
trap appears to be slightly weaker. We attribute this observa- 
tion to deformations of the laser field at such close distances 
to the large sphere. Nevertheless, the overall agreement is 
excellent and validates the theoretical approach on a quanti- 
tative level of a few percent. In particular, we conclude that 
even at the closest distance the point-particle limit is accu- 
rate not only for the stationary diffusion coefficient but also 
for the time-dependent motion. 




FIG. 6: (Color online) Normalized time-dependent diffusion coef- 
ficients, {t)/D and D±{t)/D, with the bulk diffusion constant 
D = ksT /&-Kr]a for the direction parallel and perpendicular to 
the wall at the same distances h as in Fig. |5] The experimental 
data are represented by symbols and the full lines correspond to the 
theoretical fits. 



In the experimental set-up the planar wall has been re- 
alized by inserting a much larger spherical particle in the 
fluid chamber Although theoretical models have not ad- 
dressed the Brownian motion close to a curved interface, 
some back-on-the envelope calculations may be made to 
confirm the validity of the approximation. First, the zero- 
frequency mobility or steady diffusion constant will be mod- 
ified and the suppression of mobility predicted by Lorentz, 
Eq. (|9]), should include also terms h/R, where R denotes the 
radius of the large immobilized sphere. For the experiments 
far away from the wall, h = 9.8|jm, the suppression is less 
than 20%, and h/R ~ 1/10 suggests that the curvature ef- 
fect adds another 2%. Close to the wall, h = 4.8|jm, the sup- 
pression is already 35%, yet for the curvature an additional 
1 % should be anticipated. One may also ask at what time 
scale Tcurv the effects of curvature should be manifest in the 
correlation functions. Assuming that curvature is relevant 
if the fictitious flat wall is separated from the immobilized 
sphere by more than h, simple geometric considerations lead 
to Tcurv = 2piRh/rj. In our experiments these times are at 
least a factor of 20 larger than the corresponding t„ and the 
signal is indistinguishable from noise. 



11 



VII. CONCLUSION 

Optical trapping interferometry allows monitoring the 
motion of a single colloidal particle on time scales where 
its momentum plays an important role. Due to the hydrody- 
namic memory of the fluid friction, the velocity autocorre- 
lation function (VACF) exhibits an algebraic long-time de- 
cay rather than an exponential relaxation. If the colloid is 
placed in close vicinity to a bounding wall the vortex diffu- 
sion is significantly hindered leading to a more rapid decay 
of the VACF. Although weak optical trapping is employed, 
where the trap relaxation time scale exceeds the ones of the 
fluid by orders of magnitude, the influence of the confining 
harmonic potential becomes significant. At times where the 
algebraic decay due to the wall should be visible the trap in- 
troduces additional features in the signal which have to be 
disentangled carefully. To analyze experimental data the full 
frequency dependence of the admittance including the trap 
has to be used to obtain a consistent interpretation. 

Our study should be useful to analyze experiments of col- 
loidal particles close to an interface and/or in visco-elastic 
fluids in a straightforward way. Employing the time-domain 
rather than the frequency-domain should make the interpre- 
tation of data on visco-elastic solutions, e.g. conducted by 
Atakhorrami et al |[38||, more transparent and significant, 
and improve the trapping force calibration method based on 



the power spectral density suggested by Fischer and Berg- 
S0rensen 113911 . In visco-elastic media the viscosity becomes 
itself frequency-dependent which is readily incorporated in 
the admittances. The velocity autocorrelation functions can 
be calculated numerically once the frequency dependence is 
known. Provided the experimental data exhibit little noise 
one may determine the frequency-dependent elastic moduli 
by adjusting the numerically generated curves to the experi- 
ment. Similarly, the admittances change as the properties of 
the interface change by adding surfactants, by capillary fluc- 
tuations, surface viscosity, etc. If the admittances for each 
case are known one may use the fluctuating bead as a probe 
for the local environment and determine material properties 
on scales ranging from nano- to micrometers. 
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